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Abstract 

A new method for computation of gas cooling for Lagrange approach is suggested. 
The method is based on precalculation of cooling law for known cooling function. 
Unlike implicit methods, this method is very efficient, it is an one-step method 
which is even more accurate than implicit methods of the same order. 
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1 Introduction 



The Lagrange approach in gas dynamics assumes solution of the energy con- 
servation equation in the form 

Vv--, (1) 

at p p 



where u is the internal energy per unit mass, P is the pressure, v is the gas 
velocity, A is the cooling rate per unit volume that depends on gas density p 
and temperature T. To solve (1) the implicit methods are usually used, because 
the characteristic cooling time can be much smaller than gas-dynamical time. 
This circumstance makes the usage of explicit methods impossible. 

An efficient one-step method for cooling computation for the energy conser- 
vation law in the form (1) is suggested in this paper. It is implied for the use 
in the computational gas dynamics and it is intended to predict variation of 
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the gas temperature over a time step. This method utihzes the information on 
given coohng function A(p, T) and uses approximation in which coohng time 
is shorter than gas-dynamical time. The last condition is fulfilled in many 
gas-dynamical problems, particularly, in astrophysical problems which are in 
the area of the author's interest. 



2 One-step method 



Since cooling function depends on the temperature, let rewrite the equation 
(1) in the temperature and density terms. Assuming the equation of state of 
an ideal gas with mean molecular weight /i and adiabatic exponent 7, equation 
(1) gives 



dT 



(Vv)r + f A 

Kb P 



(2) 



where is the Boltzmann constant. The first term in the right-hand side 
of this equation conforms to purely adiabatic temperature change while the 
second one corresponds to coohng process. Following Sutherland & Dopita 
[1] and Nakasato et al. [2] we suppose coohng function to be proportional to 
the squared number density A = n^A*, while A* depends on the logarithmic 
temperature only: 



dlogT 7-1 



dt 



In 10 



Vv + 



P A*(logT) 



ph 



T 



(3) 



where the definition p = fin was used. 



Despite of condition imposed on characteristic times, adiabatic cooling or heat- 
ing may contribute significantly to temperature change. To account for this 
possibility and, simultaneously, to simplify the solving procedure, let to split 
the right-hand side of equation (3) into two independent processes, adiabatic 
and cooling ones. Since cooling time is short compared to gas-dynamical time, 
the velocity divergence and density may be taken as a constants. The resulting 
approximate expression for logarithmic temperature variation is 



log^ = -T-777(Vv)i + (7, 



where a satisfies the equation 

da _ 
dt ~ 



In 10 



7-1 p 



(4) 



In 10 pk^ 



10-'^A*((7) 



(5) 



with the initial condition (T{t — Qi) — uq = logTo, where Tq is the initial 
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Fig. 1. The way to get the logarithmic temperature variation for given initial loga- 
rithmic temperature ctq and time interval t\, — fa- The predicted value of a used in 
(4) for temperature prediction. 

temperature. Solution of this equation can be written implicitly in the form 

da; 7 — 1 p 



ao 10-^A*(a;) In 10 i^ks 



t. (6) 



Equation (5) is an autonomous one. It means that the variation of a depends 
on variation of t only. It's possible to use this property of equation (5) to 
precalculate the solution of (5). First, assume the coohng function A* to be 
defined in the temperature interval [Tmin,Tmax] or, in our case, in the corre- 
sponding interval [(Jmin, Cmax]- Second, scale the time variable as 

_ 7- 1 P . /^x 
^"lnlO//A;B 

And third, bin the values of integral (6) with lower bound fixed to (Tmm and 
upper bound running from cXmin to (7„^ax' 

da; 

-Ti . (8) 



10-^A*(a;) 



This binned function will represent the cooling process over entire range of 
temperatures where cooling function is defined. 

An example of how to find a for given initial logarithmic temperature (Tq and 
time interval is shown in fig. (1). Here the cooling function of Sutherland & 

Dopita [1] (for T > 10^ K) and Nakasato et al. [2] (for T < 10^ K) was used 
(the function plotted corresponds to heavy elements abundance five orders 
of magnitude lower than the solar one). For initial logarithmic temperature 
(To the corresponding initial time t^, in binned function (8) is searched first. 
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Next, using the binned inverse fTinction, the final value of a is obtained for 
time moment, shifted to prediction time = tg, + t. This algorithm needs two 
arrays, for logarithmic temperature and for scaled time. The actual values of 
both functions can be obtained by interpolation procedure defined for these 
arrays, i.e. just two array lookups are needed to get the temperature value by 
the expression (4). 

It is interesting to estimate the error of the one-step method compared to the 
imphcit method for solution of (1) or (2). Let denote Tqs the solution derived 
by one-step method and Tim the solution derived by implicit method 



Tim — Tq 

t 



= -(7-1) 



(Vv)Tim + -^A*(logTim) 



(9) 



It's easy to show that the errors of both methods arise in the second order. 
The second derivatives of deviation of approximate solution from the exact 
value are: 



— (lnT„,-lnr) 



dlog(A7T) 



To 



TgdTcooi dlogT 



(10) 



To 



d2 
d^ 



2 (InTim-lnT) 



To 



1^1' 
^ ^gd ^cool J 



+ 



dlog(A7T) 



rLi dlogT 



+ Co 



To 



where the gas-dynamical time and cooling time are introduced: 



— = (7-l)|Vv|, 



1 P AllogTo) 



''"cool 



;i2) 



It is clearly seen that the error of implicit method is systematically larger while 
the one-step method gives exact solution in pure adiabatic or pure cooling 
limit. 



3 Conclusions 



The presented method of calculation of cooling in Lagrange computational 
gas dynamics is faster and more precise than the implicit method of the same 
order. 
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